Gravity tests in the solar system and the Pioneer anomaly 
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We build up a new phenomenological framework associated with a minimal generalization of 
Einsteinian gravitation theory. When linearity, stationarity and isotropy are assumed, tests in the 
solar system are characterized by two potentials which generalize respectively the Newton potential 
and the parameter 7 of parametrized post-Newtonian formalism. The new framework has the 
capability to account for the Pioneer anomaly while preserving the compatibility of other gravity 
tests with general relativity. 
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A number of tests of gravity have now been performed 
in the solar system and they put severe constraints on de- 
viations from general relativity. Q] However astrophysical 
and cosmological observations show anomalies, notably 
in the rotation curves of galaxies and in the relation be- 
tween redshifts and luminosities for supernovae. Since 
gravity tests agree with general relativity, these anoma- 
lies are commonly accounted for by introducing dark 
components to the content of the Universe. As long as 
these dark components are not detected by other means, 
the anomalies can also be ascribed to modifications of 
standard gravity at galactic or cosmic scales. Obvi- 
ously, any modification of this kind also has to match the 
gravity tests in the solar system. The Pioneer anomaly 
may be a central piece of information in this debate by 
pointing at some anomalous behaviour of gravity at scales 
of the order of the size of the solar system. 

The anomaly is recorded on radio tracking data from 
the Pioneer 10/11 probes during their travel to the bor- 
ders of the solar system. Q The Doppler residuals, that is 
the differences between observed Doppler velocities and 
modelled ones, vary linearly with time, for heliocentric 
distances ranging from 20 to 70 astronomical units (AU). 
Equivalently, the anomaly can be described as a roughly 
constant acceleration ap ~ 8 x 10~ 10 ms -2 directed to- 
wards the Sun. The effect has not been explained to date 
though a number of mechanisms have been considered to 
this aim, rangingfrom systematic effects to new theoret- 
ical approaches. [j, 0, The potential importance of the 
Pioneer anomaly for fundamental physics and space navi- 
gation justifies it to be submitted to further experimental 
and theoretical scrutiny. 

In the present letter, we will focus the attention on 
the key question of compatibility of the Pioneer anomaly 
with other gravity tests. The anomaly cannot be ex- 
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plained simply from a long-range modification of the 
Newton potential. If the anomalous acceleration ap is 
ascribed to such a deviation, its value is indeed too large 
to remain unnoticed on planetary tests. We show in 
the following that this problem is cured by considering 
a natural extension of general relativity. This extension 
preserves the geometrical foundations of the latter, that 
is to say the metric character of the theory, its gauge in- 
variance, the law of geodesic motion and, therefore, the 
principle of equivalence. Only the dynamical equation of 
the metric determined by the relation between curvature 
and stress tensors is modified. The motivations of this 
modification will be discussed as well as its potential phe- 
nomenological consequences in the outer solar system. 

To this aim, we may use the assumptions of linearity, 
stationarity and isotropy. As a matter of fact, the effects 
associated with non linearity of general relativity, with 
rotation and non sphericity of the Sun, are small in the 
outer solar system. More precisely, we will consider that 
these small effects are properly taken into account in the 
standard description (see for example the general rela- 
tivistic calculation of Doppler data for Pioneer probes 
and that they have the same linear effect in the extended 
theory. As a consequence, they should disappear when 
subtracting the standard result from the modified one, 
so that the potential anomalies can be calculated with 
the simplifications of a linearized theory of gravity and a 
stationary and isotropic source. 

The phenomenological consequences of the extension 
will be characterized by two potentials accomodating the 
phenomena usually associated with a long-range mod- 
ification of the Newton potential^ and an Eddington 
parameter 7 differing from unity in the "parametrized 
post-Newtonian" (PPN) formalism. [l| We will show that 
this opens free space to account for the Pioneer anomaly 
while preserving compatibility with planetary tests. The 
extended framework also has consequences for deflection 
experiments since it leads to a range dependence of the 
Eddington parameter 7. This possibility, briefly dis- 
cussed in the end of the letter, will be analyzed in greater 
detail elsewhere^ and shown to open the way to a com- 
parison of Pioneer observations and data of most recent 
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deflection experiments. lQj| 

As a first step, we now present the motivations for 
the extended framework. As already stated, the ba- 
sic geometric features of general relativity are left un- 
changed. Motions are defined as Riemann geodesies as- 
sociated with a metric tensor g py = rj^ + h^ v with r)^ v 
the Minkowski metric (signature +1, — f, — f, — f) and hp V 
a small perturbation. The Riemann tensor Rp pua , the 
Ricci tensor Rp V — r\ p ° ' Rp pva and the scalar curvature 
R = rf lv Rp V are expressed at first order in hp V which is 
conveniently written in the Fourier domain, with k the 
wavevector. The Einstein tensor Ep V — Rp V — \r\p V R is 
transverse (k v Ep V = 0) as a consequence of Bianchi iden- 
tity, and this is also the case for the stress tensor Tp V due 
to energy-momentum conservation. 

In Einsteinian theory, these two tensors are merely pro- 
portional to each other Ep V — ^G^Tp V with Gjy the 
Newton constant. But it is easy to write a more general 
linear relation between the two tensors 
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response functions cannot be fully calculated from the 
first principles. In the sequel of this letter, we adopt a 
phenomenological point of view by considering the long 
range modifications allowed by equations Q and com- 
paring their consequences to observations in the outer 
solar system. 

We now go one step further by describing the Sun as a 
point source T pa (x) — rip^ri^M^S (x) at rest in the cen- 
ter of the solar system; M is the mass of the Sun and 5 (x) 
a 3-dimensional Dirac distribution bearing on the space 
coordinates x. As a consequence of stationarity, there is 
no time variation and the frequency ko remains null in 
the following. The Einstein tensor is written in Fourier 
space in terms of the spatial part k of the wavevector 
Ep V [k] = Xp,v 0t) [k]Mc 2 . As a consequence of isotropy, 
the linear expression of E^ takes the form 
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Xiiv pa describes a momentum-dependent linear response 
of spacetime curvature to stress tensors. This relation 
preserves transversality of E pv which can be read as a 
constraint k v Xpv pa = on the linear susceptibility. In 
spite of this constraint, the susceptibility can still take 
different forms as proven by introducing projections over 
traceless and traced components ^lj 
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The factors in front of E prT in the righthand side of these 
equations are orthogonal projectors on the two sectors 
explored by transverse tensors. The component (0) has 
a null trace and it is related to the conformally invari- 
ant Weyl curvature tensor whereas the component (1) is 
related to the scalar curvature. 

Linear response functions x^ f " J are naturally pro- 
duced by radiative corrections to general relativity. (11] 
Such corrections may for example give rise to modifi- 
cations of the lagrangian which behave, at lowest or- 
der, as quadratic forms of the curvature tensors and 
which differ in the two sectors of traceless and traced 
perturbations. [l2^ In particular, electromagnetic correc- 
tions only contribute to the sector (0) as a consequence 
of conformal invariance. Effects of such corrections in 
the inner solar system have been looked for and not 
found. 0| Long range effects may also be expected, for 
example from the Sakharov argument derivinggravity 
from a kind of elasticity of quantum vacuum. |14j Any- 
way, Einstein theory is not renormalizable^^ so that the 



This constitutes a twofold generalization of Einstein 
equation which is recovered for [k] = G^ [k] = Gn- 
First, the scalar functions G^ [k] and [k] are mo- 
mentum dependent thus having the status of running 
coupling constants. | lfil| Then, these functions differ in the 
two sectors, which will turn out to be the key point for 
accomodating a Pioneer-like anomaly. 

In order to write the solution of equations Q which is 
stationary and isotropic, we use the PPN gaugeQ 



h 00 (r) = 2<S> N (r) 
h jk (r) = 2 ($jv (r) 
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The two potentials <&at and <&p depend on the distance 
r to the Sun. The first one $at (r) represents a Newton 
potential with long-range modifications. Comparison 
with planetary observations will constrain <5>at to remain 
close to its standard expression cj) = — G ^' 1 ■ The second 
potential <&p (r) promotes the PPN parameter 7 from 
the status of a constant to that of a function. As shown 
below, the presence of $p is sensed by light, so that 
deflection experiments 0] lead to constraints on $p in 
the vicinity of the Sun. We will also see that the pres- 
ence of <i>p in the outer solar system may affect eccentric 
motions, thus allowing for a possible gravitational inter- 
pretation of the Pioneer anomaly. 

At this point, it is worth adressing briefly the question 
of nonlinearity of gravitation theory accounted for by a 
parameter (3 in the PPN metric, [lj Its effect is known 
to be significant in the inner part of the solar system 
by affecting the perihelion precession of Mercury as well 
as the polarization of the Moon/Earth orbit by the Sun. 
Its discussion requires a non linear treatment of gravita- 
tional perturbations, which is clearly out of the scope of 
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this letter. The present discussion will thus be restricted 
to the category of phenomena correctly accounted for in 
the linearized treatment. This category corresponds to 
probes having kinetic energies v 2 /c 2 much larger than 
their gravitational energies \<f>\, thus including the case of 
Pionccr-like probes escaping the solar system as well as 
the case of light. 9] 

The two potentials &n,p are directly related to the 
running gravitational constants. Writing down equations 
@ in the PPN gauge, we indeed obtain 



A$q (x) = G a (x) M 
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, a = N,P 
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where A is the spatial Laplacian operator and Gn,p are 
determined by G^ ' 1 ) 
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The standard Poisson equation is recovered when Gn — 
G N and G P = (ie = G« = G N ). In the general 
case, the functions Gn,p are momentum dependent. 

In a simple version of the framework, the two poten- 
tials contain contributions linear in r besides the ordinary 
contributions scaling as ^ 

.. (r) __2£ + fe«r , a = N , P (7) 

rc c 

This corresponds to the following running constants 
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Gn is identified as the effective Newton constant in the 
inner solar system while the 3 parameters Gp, Cn and 
£p measure the deviation from general relativity; (jv de- 
scribes a long range modification of Newton law which 
has to remain small to fit planetary data; Gp corresponds 
to the anomaly (7 — 1) of the PPN formalism (see below); 
Cp is the newest feature of the extended framework, lead- 
ing in the following to a Pioneer-like anomaly. 

It is worth noticing that expressions (|7I8|I cannot be 
exact for all distances or all momenta. What is needed 
for our purpose is that they are correct descriptions of 
the true running constants at the scales involved in the 
tests. In particular, the linear increase of J7J) with r has 
to hold roughly from 20 to 70 AU in order to explain 
the Pioneer anomaly (see below). But the same increase 
continued at much larger scales would have unwanted 
consequences for galactic astrophysics. This implies that 
expressions (J7J) can only be approximations in the solar 
system of functions better behaved at r — > 00. For exam- 
ple, they can be considered as an expansion at r <C r\ of 
Yukawa functions with a range r\ larger than the size of 
the solar system. Accordingly, the infrared divergence of 



the running constants (JSJ) is cured by an infrared cutoff 
at a wavevector of the order of l/r\. This regulariza- 
tion does not solve the problem of the matching of solar 
system physics with larger scales but it is sufficient to 
discuss phenomena within the solar system. Anyway, (J7J 
is only a simple model for the potentials $jv,p with the 
advantage that the phenomenology is determined by a 
small number of constants but the drawback of a loss of 
generality with respect to the more general r-dependence 
allowed for the metric Q. 

We now discuss the potentially observable conse- 
quences of this framework. To this aim, we study 
geodesic motion in the metric (@J or, equivalcntly, 
Hamilton- Jacobi equation for wave propagation. We de- 
note by t the time coordinate, r the radius, ip the az- 
imuth angle, 9 the colatitude and suppose the trajectory 
to take place in the plane 9 = \. For matter waves 
with a non null mass for example, we use the conserva- 
tion of energy E = mc 2 j o^ and angular momentum 
J = mr 2 sm 2 9gn^ where ds is the invariant length el- 
ement along the motion. We also denote by v r = 



and 



dip 

ds 



the velocities measured respectively in 



the radial and orthoradial directions. In a first stage, we 
briefly discuss the effect of a modification of the first po- 
tential <E>jv, setting <E>p to zero, and we recover the known 
fact that it cannot account for the Pioneer anomaly. We 
then shift our attention to the effect of $p alone. In the 
concluding paragraphs, we will sketch the program of re- 
analyzing gravity tests in the general case where $ n can 
be modified and <E>p differ from zero. 

A long range modification of the Newton potential $ n 
could be detected as an anomaly of the third Kepler law 
on a circular orbit (v r — 0). To evaluate this effect, we 
compute the square l> 2 of the orthoradial velocity in the 
modified theory and subtract from it the constant value 
[w 2 ] st obtained in standard theory. We denote by 5v 2 
the difference which measures the potential anomaly 



5v* = v\ 



VJst 



( N Mr 



(9) 



This anomaly, which is proportional to Cat, is not ob- 
served on the motions of planets or probes in the solar 
system. In particular, the telemetry data on probes 
close to Mars are sufficient to put an upper bound on 
Cat which is much too small to account for the Pioneer 
anomaly (see section XI-B in Q). After a translation 
into the notations of the present paper, the bound on 
the anomalous acceleration |C/vM| < 5 x 10~ 13 ms~ 2 is 
found to be much smaller than ap. Since a long-range 
modification of the Newton potential $ n cannot explain 
the anomaly, we now consider the case where the first po- 
tential $tv has its standard form and focus our attention 
of the effects of the second potential $p. 

The potential $p affects Dopplcr tracking data of 
Pioneer-like probes with eccentric motions. To evaluate 
this effect, we calculate the motion of such probes in the 
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metric (@J and also take into account the perturbation of 
the propagation of radio signals to and from the probes. 
We then express the result as an equivalent acceleration 
a defined as the time derivative of the Doppler velocity. 
We finally obtain the anomaly by subtracting the stan- 
dard expression calculated with general relativity. The 
anomaly is found to be proportional to the derivative 
of the second potential as well as to the square u 2 
of the radial velocity of the probes^ 

d$ P , / G P M\ vl , s 

6a ee a-[a] st ~ 2^z, 2 ~ 2 ^ P M + -^-j -| (10) 

We know that Gp <C Gn (see also below) and u 2 <C c 2 , 
so that the term proportional to Gp can be neglected 
in (|10|l . We are therefore left with a constant anoma- 
lous acceleration Sa directed towards the Sun, if £p has a 
positive sign. Tentatively identifying this result with the 
Pioneer anomalous acceleration fixes the unknown pa- 
rameter \C,pM\ = a P c 2 /(2vf) ~ 0.25ms~ 2 which is thus 
found to be much larger than |£jvM|. 

Equation (|10|l means that the new framework pre- 
sented in this letter effectively has the capability of ac- 
comodating a Pioneer-like anomaly for probes having a 
large radial velocity. It also leads to the prediction that 
the anomalous acceleration should show a dependence 
versus the velocity of the probes. As the two Pioneer 
probes have nearly equal velocities and nearly equal ac- 
celerations, this prediction cannot be confronted to avail- 
able data. But it could be checked out by confronting 
(|TU|l to the data recorded on Pioneer probes. This pre- 
diction also has to be kept in mind when proposing new 
missions, since it points to the idea of trying probes with 
different radial velocities. Equation (|10l) also predicts 
a specific r— dependence for the anomalous acceleration 
when the term proportional to Gp and the variation of 
u 2 on the trajectory are kept. This might be tested not 
only in the outer solar system but also on probes flying 
to Mars or Jupiter, if the sensitivity of the acceleration 
measurement can be made good enough in spite of per- 
turbations such as solar wind and radiation pressure. 
This might be done by embarking an accelerometer in 
the probe and thus distinguishing between drag forces 
and geodesic motion. 

In the sequel of the letter, we discuss other phenomeno- 
logical consequences of the presence of $p and check out 
that the modification of Einstein theory needed to obtain 
(|10|l does not spoil its good agreement with other gravity 
tests. A critical problem in this context is the effect of 
$p on propagation of light rays. In order to assess this 
effect, we compute the deflection angle tp for light rays 
passing near the Sun in the metric Q and subtract the 
standard value to obtain the potential anomaly as0 

m - i r/i 2GpM 2 ^ pM P r , s 

5V = V " Met -2 ^ L(P) 



p is the impact parameter, that is also the distance of 
closest approach of the light ray to the Sun; L is a factor 
of order unity which depends logarithmically on p and 
on the distances of the emitter and receiver to the Sun; 
for an emitter far outside the solar system, its distance is 
replaced by the range r\ at which the linear dependence 
of the metric falls down to zero. Should £p be set to 
zero, equation lfTT)l would be equivalent to the PPN re- 
sult (see 0) with an Eddington parameter 7 determined 
from the ratio of gravity constants Gp/Gn = (1 — 7). 
Deflection tests would thus tell us that Gp/Gn is much 
smaller than unity with a maximum value given by the 
upper bound on (7 — 1). The novelty in equation 
is the term proportional to £p which implies that de- 
flection tests could show an anomaly depending on the 
impact parameter. This result certainly pleads for a re- 
analysis of the most recent deflection data obtained dur- 
ing the Cassini mission. [T^j A comparison of these data 
with or, preferably, the more precise expression in 
0, would indeed lead to a test of the presence of Cp in- 
dependent of the Pioneer test. An observation of a non 
null (p would constitute a clear indication for the new 
framework. Note that different values of £p could also be 
deduced from Pioneer and Cassini data since Cp is in fact 
a function which could vary when going from the radius 
of the Sun ~ 0.7 x 10 9 m to the distance explored by the 
Pioneer probes ~ 1.2 x 10 13 m. On a longer term, the /in- 
dependence of deflection experiments also constitutes a 
further motivation for high accuracy deflection tests such 
as LATOR 17j or astrometric surveys such as GAIA.fljf 
The present letter constitutes a first study of the phe- 
nomenological consequences of the extended framework. 
As already stated, the modified gravity equation natu- 
rally leads to metric perturbations characterized by two 
potentials <&at and $p. It is therefore necessary to per- 
form a new analysis of the motions in the solar system 
looking for the combined effects of these two potentials. 
The nonlinearity of gravitation theory has also to be 
taken into account in order to evaluate for example the 
precession of perihelion of planets. Clearly the eccentric- 
ity of the orbits will play a key role in these discussions; it 
takes large values for Pioneer-like probes which sense <J> p 
whereas it is zero for circular orbits which do not. This 
suggests to devote a detailed analysis to the intermediate 
situation, not only for the two categories of bound and 
unbound orbits, but also for the flybies used to bring 
Pioneer-like probes from the former category to the lat- 
ter one. 0,0] 
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